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Based on renormalization group arguments we establish that for a superconductor in the pres- 
ence of a weak external magnetic field, B, the dependence on B and the deviation from the critical 
temperature, r, of a thermodynamic quantity, P, takes the scaling form P = t e X {-§^r~ 2v , qr~ vu ' e ), 
where 9 and v are XY exponents, q is the scaled electromagnetic coupling and vu^ is the associated 
crossover exponent. For g/r"" c <C 1, the experimentally accessible region in high-T c superconduc- 
tors, there is a reduction to one-variable scaling plus small corrections. In this region we find the 
shift in the specific heat maximum is given by A = xo(B /Qq) 1 ^ 2 " and that the singular part of 
the free energy at the critical temperature takes the form F a i ng — c(d)(_B/$o) d ^ 2 where c(d) is a 
universal amplitude. A one loop approximation in three dimensions gives c(3) ~ 0.22. The results 
presented here should have equal applicability to the nematic to smectic-A transition. 
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Since the discovery of high-Tc superconductivity there 
has been a significant increase in interest in the critical 
properties of superconductors. As one of the principal 
experimental probes for the analysis of high-T c materi- 
als involves their study in external magnetic fields it is 
important to have reliable and robust theoretical predic- 
tions for such a setting. 

The formal analogy between a superconductor in an 
external magnetic field and a system constrained by fi- 
nite size is well known ffl. In Q we established that 
correlation functions of such a superconductor exhibit 
a scaling form identical to that of a near critical sys- 
tem in a confined geometry of size L, with the combi- 
nation (£?/$o)~ 1//2 , where B is the magnetic field and 
<&o = itTicje is the unit flux quantum, playing the role 
of L. A crucial ingredient in our analysis was the as- 
sumption of a non-trivial, infrared stable renormalization 
group (RG) fixed point. In conventional superconductors 
such a fixed point is experimentally inaccessible. Ad- 
ditionally, there are theoretical arguments || that cast 
doubt on the existence of such a fixed point except in the 
artificial case where the number of components of the or- 
der parameter is large. For type-II superconductors there 
is some evidence Q that very close to T c the relevance 
of magnetic field fluctuations causes a crossover to a sta- 
ble "charged- XY" fixed point with inverted temperature 
axis. Once again in the case of conventional supercon- 



ductors this crossover would be unobservable. 

On a more pragmatic note, experimentally there is now 
ample evidence ||-|7j that for high-T c superconductors, 
for |T — T c \ < 10K, the zero field transition exhibits 
critical behaviour characteristic of a three dimensional 
AY-modcl. This fact, together with the finite size anal- 
ogy, has been utilized Q to collapse experimental spe- 
cific heat data using a "finite size" scaling ansatz of the 
form P = t 6 X{J^t~ 2u ), where r is a linear deviation 
from the critical temperature and v and 6 are XY expo- 
nents, with v ~ 0.67 for a three dimensional sample. In 
a recent letter Lawrie Q presented two-loop arguments 
establishing that to this order, when fluctuations in the 
electromagnetic field are neglected, the two-point func- 
tion of the order parameter in a magnetic field exhibits 
such a scaling form. 

The neutral- XY fixed point where e* = is unsta- 
ble with respect to a small perturbation of the electro- 
magnetic coupling. Hence, one expects a region wherein 
electromagnetic fluctuations can be treated as a relevant 
perturbation, in a similar fashion to temperature, thus 
leading to a two-variable scaling formulation. For high- 
Tc superconductors, where the Ginzburg-Landau param- 
eter, k, is large (e.g. for YBa2Cu 3 07- y , k ~ 87 the 
corrections due to electromagnetic fluctuations, although 
relevant, are very small except in a region asymptotically 
close to T c and can be taken into account as corrections 
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to scaling around the single variable scaling form asso- 
ciated with the neutral-XY" fixed point. Experimental 
evidence certainly seems to accord with this picture. 

In this letter we establish using RG techniques that 
the appropriate scaling form for a general thermody- 
namic quantity P, in a magnetic field takes the form 
P = t e X{-§^T- 2v ,qT~ VWe ). In the case where electro- 
magnetic fluctuations are weak, such as high-T c super- 
conductors, we show that this reduces to a one-variable 
scaling form plus small corrections to scaling relative 
to the XY fixed point. In addition, our analysis pro- 
vides several new theoretical predictions which should 
be testable as experiments on high-T c superconductors 
become more refined in the vicinity of the critical point. 
The derivation elucidates the assumptions on which such 
a scaling form relies, in particular the result is quite ro- 
bust. 

Our starting point is the Landau-Ginzburg- Wilson 
Hamiltonian for a superconductor, based on the Cooper 
pair density as order parameter, given by 

H[<p] = /d d x(j(V - ff AM 2 + roM 2 + f M 4 

+ 1(VxA) 2 -IB.Vxa), (1) 

Note that with the substitutions, r = A = ^j^- 
and ip — ip we recast Eq. ([j]) into conventional form 
(as found in e.g. Blatter et al. JlOfl). The only tempera- 
ture dependent parameter is a (and hence ro), which is 
assumed to depend linearly on T . Though high-Tc su- 
perconductors are highly anisotropic, if the anisotropy 
is not so strong as to invalidate a continuous LGW for- 
mulation, the asymmetric results can be obtained from 
the symmetric ones by the rescaling methods of Blatter, 
Gershkenbein and Larkin pi| , hence we concentrate on 
the symmetric case. 

In the case of a superconductor in a unidirectional 
magnetic field, B, fluctuations of ip induce a renormal- 
ization of B even if electromagnetic fluctuations are ig- 
nored. This renormalization becomes ultraviolet diver- 
gent at the upper critical dimension and is encoded in 

the wavefunction renormalization of the vector potential 
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A, i.e. A = Zjl A , (where we use a subscript to 
indicate a bare or unrenormalized quantity). The elec- 
tromagnetic coupling e = Z e e$ is similarly renormalized. 
Gauge invariance as implemented by the Ward identities 
relates the two renormalizations, ZeZ 1 / 2 — 1 fH| . An 
immediate consequence is that eB is an RG invariant. 
Furthermore, from the Hamiltonian (Q) we see that only 
the combination P/$o arises in any correlation function 
or thermodynamic quantity associated with the ip field. 
Thus {Qq/B) 1 / 2 is a nonlinear scaling field with the di- 
mensions of length which behaves for RG purposes in a 
fashion very similar to the linear size I of a finite or 
layered system. 



We begin by establishing the scaling form for a gen- 
eral correlation function. From the Hamiltonian ([!]) we 
see that correlation functions depend on ro, Ao, bq and 
the RG invariant B/Qq together with a cutoff A. Defin- 
ing the renormalized correlation functions G^ N,M ' by 

< p*(l)..p*(k)p(k + l)..p{N)\p(l)\ 2 ..\p{M)\ 2 > 

= G {NM) (t, A,e, — ,fj,,xi, .,x N ,yi, ;Vm) (2) 
<±>o 

where t is the renormalized deviation from the zero field 
critical temperature, i.e. ro = r c + Z^it, the relation 
between the bare and renormalized correlation functions 
is 

G™(t, A, e, M ) = Z^Z^ N - M \r , A , e , A) 
where for convenience we have suppressed the depen- 
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dence on position. Z v and Z^i are the field and com- 
posite operator renormalization factors, Z\\ = Ao and fi 
is an arbitrary renormalization scale. Here we take A and 
e to be the dimensionless renormalized couplings. 

The Hamiltonian ([!]) defines a renormalizable field the- 
ory for any dimension d < 4. As ultraviolet divergences 
are unaffected by the presence of an external B field the 
Zi obtained via a B independent renormalization scheme, 
such as minimal subtraction, are sufficient to renormal- 
ize the theory for B ^ 0. The RG equation obeyed by 
q(n,m) j s e q U i va i en t to the statement that the bare func- 
tions are independent of the arbitrary renormalization 
scale chosen to define the renormalized theory; hence the 
renormalized correlation functions obey the equation 

~(M V + ^7„)C<«« = 0, (3) 

where f3 e = fide/dfi, (3\ = /idX/dfj, and 7 V = ^idliiZ v /dfj, 
and 7^2 = —fidlnZ^-i/dfi are the anomalous dimensions 
of (p and (p 2 respectively. The solution of Eq. (|J) is 

G (M) Mi 4 irf=e /;(^w+«vw)? 

xG^ M \t(p),X(p),e(p),^-,pfj,). (4) 

The absence of a term of the form psd/dB in (||) is due 
to the fact that B only appears in the RG invariant com- 
bination i?/$o and is directly analogous to the absence 
of a term Phd/dL in the RG equation for a finite system. 
The latter is a crucial element in the proof of finite size 
scaling Jl3| . 

To proceed further we need to make some assump- 
tion about the fixed points of the model. The Landau- 
Ginzburg model under discussion has several potential 
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candidates. In @ we assumed an infrared stable fixed 
point, e = e*, A = A*. By going to a "dual" formulation 
there is some theoretical evidence j| that for type-II su- 
perconductors the zero field transition is associated with 
a fixed point of this type which is in the universality 
class of the charged-XY-model with inverted tempera- 
ture axis. This fixed point is associated with only one 
relevant direction, t, and in its vicinity scaling functions 
depend on only one scaling variable. In the case of super- 
conductivity there is at present no experimental evidence 
for such a fixed point. Another candidate is the Gaussian 
fixed point A* = 0, e* = 0, however, a scaling formula- 
tion based on an expansion around this fixed point is in 
disagreement with experimental results for high-T c su- 
perconductors. 

Given that there is compelling experimental evidence 
that the XY fixed point is relevant to critical fluctua- 
tions in high-T c superconductors we will now concentrate 
our attention on a scaling description based on the fixed 
point: e* = 0, A = A* where A* is the Wilson-Fisher fixed 
point coupling for an 0(2) model. One might worry that 
at such a fixed point, since e* = 0, the dependence on B 
would drop out. This is not the case since eB is an RG 
invariant and if it is non-zero at the beginning of an RG 
flow it is non-zero at the end. 

There is an essential difference between this fixed point 
in the context of a superconductor and the standard XY 
case since in the scaling neighbourhood of this neutral- 
XY fixed point, as can be seen from the work of Halperin, 
Lubensky and Ma ||, there are two relevant variables, 
t and e, and therefore one has generically two- variable 
scaling functions, with t being the leading relevant per- 
turbation. As /3 e (A, e) = — w e (A, e)e linearization about 
e* = ensures that there is no mixing of e and A in 
the relevant direction, though the irrelevant direction is 
a linear combination of e and A || . For d = 3 estimates 
give v ~ 0.67 and u> e ~ 1. 

Since the RG equations for the two relevant variables 
are 



de 
dp 



-u) e (\,e)e and p^- = j v 2(e, X)t, (5) 



in the scaling neighbourhood of the fixed point we have 
7y> 2 (A*,0) = 2—1/ is, where v is the XY correlation length 
exponent and w e (A*,0) = lu e is the scaling dimension of 
the charge. Integrating (||) yields t(p) — tp 2 ~ 1 / 1 ' and 
e(p) = ep'^". The general scaling form of (||) then be- 
comes 



G' 



'N,M r>{N,M) 



B 



Q 



(6) 



where 9 NM = §(vd - 7) + M(ud - 1), r = tp 1 ^- 2 , 
q = e/i" e and G( N > M > has been rescaled by a fac- 



tor p 



-f v -M(2-i) 



Since the co-ordinates Xi do not 



g(N,M) m combination Xipp. Eq. (^J) implies G^ N,M ^ 
is a homogeneous function, therefore the variable p can 
be used to eliminate one of the arguments. We choose to 
eliminate p by setting r/(p/x) 1 ^ = 1. The general scal- 
ing form in the vicinity of the neutral-AY fixed point is 
then 



Q{N,M) =T 6 N , M g{N,M)^ y) 



(7) 



where the scaling variables are x = B/^ot 2ij and y — 
q/T VU3 ° and 0(^*0 is a universal scaling function. A 
single-variable scaling function is recovered in the limit 
q — > 0, or alternatively in the limit T — > T c if there exists 
a non-trivial stable infrared fixed point. 

The free energy takes the two- variable scaling form 



- a F{x,y) 



(8) 




with x and y as above. When the temperature is precisely 
tuned to the critical temperature it is no longer appro- 
priate to eliminate p in terms of the temperature variable 
t but instead one can use the condition B /&o(pp) 2 = 1. 
The dependence on the external field of the free energy 
then takes the simplified form 



(9) 



with a corresponding expression for the specific heat. 

We have shown here that generically one expects to 
see two-variable scaling functions when electromagnetic 
fluctuations are taken into account, however, one- variable 
scaling seems to give a very good fit to current experi- 
mental data. The consistency of these two seemingly 
contradictory facts can be seen by making a Taylor ex- 
pansion of the scaling functions (0-^)- One finds, for 
instance, for the correlation functions to first order in e 



q{n,m) _ T e N ,M(i _|_ 



qr 



'D{x))Gt N > M \x). (10) 



get renormalized they enter the dimensionless functions 



For strongly type-II superconductors the correction to 
scaling factor, even though it is associated with a relevant 
operator, will be very small unless one gets very close to 
the critical point, and hence a one- variable scaling form 
associated with the XY fixed point should offer a very 
good description. Naturally it would be very interesting 
to see if such corrections to scaling due to electromagnetic 
fluctuations could be observed experimentally. 

With the scaling formulation given above and the anal- 
ogy with finite size scaling there are several interest- 
ing corollaries regarding critical temperature shifts and 
rounding. For d < 4 one does not expect any critical 
divergence for B ^ 0, hence any shifts will be associated 
with a pseudo-critical temperature T P (B), which we take 
to be the point at which the specific heat is a maximum, 
i.e. dC/dr\ T=Tp(B) = 0. Given that the specific heat is 
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proportional to J d d xG^ ' 2 ^ 
scaling form 



then r(°' 3 > = dC/dr has the 
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r^(T,\,e,—)=T d »- 3 M(x,y). 

<J>0 



(11) 



The shift in temperature of the specific heat maximum 
A oc T P (B) — T c (0) is given by the zero of the uni- 
versal scaling function Ai, i.e. Ai(x, y) = 0. Once 
y = q/A.^' « 1 we can set q = in (|ll|), a good ap- 
proximation experimentally, (corrections to scaling can 
of course be added). Denoting this universal zero, xq, 
one finds 



l/2v 



(12) 



which is identical to the corresponding formula familiar 
in finite size scaling fl4|| . A similar expression for the 
specific heat rounding follows from the standard argu- 
ment that finite size effects will start to become impor- 
tant when £o ~ {B/^a)^ 1 ^ 2 , where £o is the zero field 
coherence length. 

The above general argument can also be used to derive 
in the present context the analog of phenomenological 
renormalization p5| . Defining the coherence length via 
the second moment of the two-point correlation function 
leads to 



Co V $0 



(13) 



where 72. is a universal scaling function and once again 
we are neglecting electromagnetic corrections to scaling. 
Defining a scale transformation B — > B' , t — > t' via 
B£ 2 (t,0) = B'f 2 (i',0) leads to 



e It, 



B 



ie a 



B'/B 



(14) 



which represents an exact RG transformation. T = T c (0) 
is a fixed point of this equation and the critical exponent 
v can be obtained by linearizing around this fixed point 
to find 



ML 

<I>o 



(15) 



In the approximation where electromagnetic fluctua- 
tions can be neglected, or if a charged-XF fixed point 
exists, the above scaling formulation also leads to many 
new universal critical-point ratios. For instance, the ra- 
tio £ 2 -B/$0: evaluated at the critical temperature, should 
be a universal number. Also, the singular part of the 
free energy at T = T c (0) will have the form F S i ng = 
c{d)(B/<f> ) d/2 where c(d) is a universal number analo- 
gous to the Casimir amplitude familiar in finite size scal- 
ing pq| . A one loop approximation gives 



= -r(i - ~)C(i - ~ 5) 



(16) 



which for d — 3 gives c(3) ~ 0.22. 

To conclude: we have derived using RG methods the 
scaling behaviour of thermodynamic functions in the crit- 
ical region for a superconductor in a magnetic field. We 
find that in the neighborhood of the XY fixed point the 
relevance of electromagnetic fluctuations leads to a two- 
variable scaling formulation. For q/T l,UJ " <C 1, the exper- 
imentally accessible region for high-T c superconductors, 
this two-variable formulation reduces to a one-variable 
form plus corrections to scaling due to the electromag- 
netic fluctuations. 

We believe that various predictions we have made here 
should be experimentally verifiable in high-T c supercon- 
ductors. In particular, the analogy with finite size scaling 
is robust and it should be possible to determine in exactly 
what regions the data collapses onto a one- variable scal- 
ing formulation and to what extent corrections to scaling 
can describe any deviations. It should also be possible to 
measure some of the universal critical point ratios that 
arise naturally in the present scaling formulation and the 
predicted shift of the specific heat maximum. It would 
be interesting to extend the present approach to include 
dynamics, anisotropy and/or disorder. The conclusions 
of this letter should also be of relevance to the nematic 
to smectic-A phase transition in liquid crystals which is 
believed to be in the same universality class. 
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